Gauge fields, ripples and wrinkles in graphene layers. 
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We analyze elastic deformations of graphene sheets which lead to effective gauge fields acting on 
the charge carriers. Corrugations in the substrate induce stresses, which, in turn, can give rise to 
mechanical instabilities and the formation of wrinkles. Similar effects may take place in suspended 
graphene samples under tension. 

PACS numbers: 73.20.-r; 73.20.Hb; 73.23.-b; 73.43.-f 



Introduction. Two dimensional graphene membranes 
define a new class of materials with a number of novel and 
promising features P, 0, El- From a basic point of view, 
one of the striking properties of graphene is that many 
types of long wavelength disorder influence the electrons 
as effective gauge fields 0, i, H 0, 1 i, M, El E El G5, 
TBI ], which modify the electronic properties, as measured, 
for instance, in weak localization experiments [?], d, [H, 

13- 

Because of the existence of two inequivalent valleys 
in the Brillouin Zone of graphene, one can define for- 
mally two types of gauge fields, one of which induces 
mixing between the valleys, and another which does 
not. These two fields do not commute. A gauge field 
which effectively rotates electrons from one valley into 
electrons in the other valley is needed to describe lat- 
tice disclinations, such as pentagonal and heptagonal 
rings 0, H, [3, [3]. Intravalley gauge fields can be in- 
duced by elastic strains fH H2I liolTill] . by curvature [H, or 
by lattice defects such as dislocations. For each type of 
field, the hamiltonian can be written as two independent 
Dirac equations, with gauge fields of opposite sign for 
each of them, so that it is invariant with respect to time 
reversal. 

Elastic strains of reasonable values, u ~ 0.02 — 0.05 
with the correct symmetries can lead to significant ef- 
fective magnetic fields, B ~ 0.1 — 1T0, EH. Large 
graphene sheets show corrugations (2^, and significant 
scale deformations (23l 24 , 25 1 . Suspended samples [25II2II] 
can also be deflected if a gate potential is applied to 
them[23|. Graphene deposited on a substrate shows 
corrugations [3, 0, H^, The existence of corruga- 
tions is confirmed by numerical simulations [3lj| . 

The low energy electronic states of graphene are de- 
scribed by the Dirac equation^. The properties of the 
electronic spectrum of the Dirac equation in the pres- 
ence of a random gauge field has been studied in connec- 
tion with the Quantum Hall Effect [32I I33I l3~j] . The scat- 
tering cross section of isolated ripples vanishes at short 
and long electron wavelengths, with a peak for wave- 
lengths comparable to the size of the ripple [351] . The 
density of states is changed at low energies, and it ac- 



quires an anomalous power law dependence on energy, 
where the exponent depends on the strength of the dis- 
order. The combined effect of random gauge fields and 
the long range electron-electron interaction leads to a 
complex phase diagram when the system is close to the 
neutrality point 36 
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In the following, we discuss some basic features of rip- 
ples and wrinkles in single layer graphene. We analyze in 
the following section the elastic strains in a periodic array 
of ripples, similar to that reported in [3^]. and the effect 
of the induced gauge fields. We next discuss possible 
mechanisms which may induce wrinkles [40| in graphene, 
and the way they can alter the electronic properties. 

We will not consider here the effect of a scalar poten- 
tial also induced by the presence of strains [20j]. When 
the length scale of the corrugations is larger than the 
Thomas- Fermi screening length k~ x ~ kp 1 , where fcp 
is the Fermi wavevector, these scalar fields will be 
screened [27j| . while the gauge fields are not. This pic- 
ture is confirme d by full electronic calculations of rippled 
graphene layers [14j. 

Periodic array of ripples. The simplest corrugation 
which leads to inhomogeneous strains and, as a result, to 
non trivial gauge fields in graphene is a periodic array of 
ripples, which can be induced by a suitable substrate [39(. 
For simplicity, we assume the existence of modulation of 
the out of plane displacement of the graphene layer of 
the type: 

h(x, y) = h [2 cos(Gir) + 2 cos(G 2 r) + 2 cos(G 3 r)l (1) 

where G*i = Gn x ,G 2 = G(n x /2 + %/3n y /2) and G 3 = 
G(— n x /2 + v3n y /2). These vectors define a triangular 
lattice. In terms of the in-plane displacements, u(x,y), 
the strain tensor is: 



l yy ~ ®y u y 



2 

{dyhf 



d y u x + d x u y (d x h) (d y h) 



(2) 
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FIG. 1: (Color online). Top: Compression, u xx + u yy , induced 
by a periodic array of ripples, see text for details. Bottom: 
Effective magnetic field for the same array. 



Neglecting for the moment the contribution due to the 
bending stiffness of the graphene layer, the free energy 
can be written as 4l| : 



.?"=-/ d z r{u xx 



'yy 



d 2 f i 



2u 



"yy ' °""®y) 
(3) 

where, using for the l ong itudinal and transverse sound 



j, |44| the values v L = 22.2 
14.7 Km/s, we find A = 1.6eV A" 2 and 



velocities of graphene [42 
Km/s and vt = 
fi = 5.7eV A" 2 . 

The full strain tensor tiy can be obtained from the 
tensor = (dih) (djh) by minimizing the free energy, 
e q. (|3l) , and obtaining the in plane displacement field, 
m[13(. In terms of the strain tensor, the effective gauge 
field acting on the electrons, for a given valley in the 
Brillouin Zone, isj^ij: 



A X {r) — — ( U XX Uyy) 



a dlog(a) 
, c$q d\og(t) 
a d\og(a) 



(4) 



where (3 = d log (t) /dlog (a) w 2 gives the dependence 
of the tight binding hopping element between orbitals 
in nearest neighbor carbon atoms, t 3eV, on the dis- 
tance between them, a ~ 1.4A, c is a numerical constant 
of order unity, and <&o is the quantum unit of flux. In 
this analysis, we have assumed that the shape of the cor- 
rugation, h(x,y), is determined by external forces. We 
now consider the bending rigidity of the layer. From the 



FIG. 2: (Color online). Density of states for a periodic array 
of ripples, including the induced gauge field. Blue: ho = 0. 
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continuum theory of elasticity 4lJ we expect a term: 
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(5) 



For graphene, n\ ~ ki ~ R ~ leV. The ratio between 
the bending and compression energies of the ripples is 
approximately proportional to G 2 R/(X + For meso- 
scopic ripples of lengths I — 2irG~ 1 ~ 10 — lOOnm, the 
bending rigidity can be neglected. 

The calculation of the strain tensor and effective gauge 
fields for the height profile in eq. ((T|) is rather simple. Val- 
ues for the strain and induced magnetic field are shown in 
Fig- [2 The electronic density of states, calculated includ- 
ing the gauge field due to the ripples, is shown in Fig. [2] 
The results have been obtained by using a discrete lattice 
model with 2x81x81 sites with a broadening of the levels 
of Sui — 0.03t. The small amplitude undulations in the 
density of states are an artifact due to the discretization. 
There is a peak at the Dirac energy for the largest corru- 
gations, consistent with analytical results for the density 
of states of Dirac electrons in the presence of a random 
gauge field [34j. Similar effects are obtained in numerical 
studies of the Dirac equation in a random field whose 
fluctuations are smooth on the scale of the lattice jiU, al- 
though they have not be found if the gauge field varies 
on length scales comparable to the lattice spacing (33l|. 
The highest effective magnetic fields induced by ripples 
with ho ~ lnm and I — 2irG~ 1 ~ 60nm can be large, 
B m ax ~ 2 — 5T[7], so that the number of flux quanta per 
ripple exceeds one. 

Wrinkling of graphene layers: Graphene on a sub- 
strate. As in any membrane, topological lattice de- 
fects can spontaneously lead to corrugations [46| in 
graphene (see also Sec. VIB in [13|). Similar effects 
can be expected from impurities which induce lattice 
deformations Rigid layers on soft substrates show 
wrinkles if the are under anisotropic stresses which ex- 
ceed a given threshold (iol 48 . [49| . A related situation 



FIG. 3: (Color online). Generation of wrinkles in a graphene 
layer on Si02- Top: The layer is compressed against the 
corrugated substrate. Bottom: When the external pressure is 
removed, in-plane stresses due to the deformation of the layer 
induce wrinkles, if the pinning to the substrate is sufficiently 
weak. 
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FIG. 5: (Color online). Top: Wrinkles expected in a sus- 
pended graphene sheet of I = 200nm under an applied tension 
Al/l = 2 x 1(T 2 . Bottom: Effective magnetic field induced 
by the wrinkles shown in the top. 



FIG. 4: (Color online). Stresses along the line (a;, 0) for the 
periodic array of ripples analyzed in Fig. [TJ Red: a xx (x,0). 
Green: a yy (x,0). 



which may be relevant for graphene layers on SiC>2 sub- 
strates is sketched in Fig. |3] The exfoliation procedure 
applies pressure on the graphene layer against the cor- 
rugated substrate, inducing deformations and in plane 
stresses. When the pressure is removed, the system can 
reduce the elastic energy by forming wrinkles, provided 
that the stresses are sufficiently large and the pinning 
energy to the substrate is sufficiently low. 

We show in Fig. [4] the stresses a xx and a yy induced 
by the periodic array of ripples discussed in the previous 
section. There are regions, whose length is fraction of 
G _1 , where the stresses are highly anisotropic. The gain 
in elastic free energy per unit area due to the formation 
of wrinkles in these regions is T ~ (A + 2^,)(Gh ) 2 ~ 
0.05 - O.leV k~ 2 , for ho ~ lnm and I = ItxQ- 1 ~ 60nm. 
This value is larger than typical interactions energies be- 
tween a graphene layer and a SiC>2 substrate, or other 
materials, like water or ions, which may be trapped be- 
tween graphene and the substrate [Boj]. Hence, we ex- 
pect that, if graphene is pressed onto a corrugated SiC>2 
substrate, the removal of that pressure will lead to the 
spontaneous formation of wrinkles. Note that the low- 
ering of strains will actually reduce the effective gauge 
field acting on the electrons. On the other hand, the 



regions detached from the substrate can support low en- 
ergy flexural modes, which can scatter electrons at finite 
temperatures (Hli [o2|. 

Wrinkling of graphene layers: suspended graphene 
samples. We now analyze the structural instabilities to- 
wards curved shapes, wrinkles, which can arise in sus- 
pended membranes under tension[53j, see also Sec. Vll 
in [13| . A stretched membrane with clamped ends will de- 
velop wrinkles if the applied tension is large enough (40j. 
This threshold arises from a balance between the bending 
energy and the applied tension Wrinkles will occur 
when the tension, T ~ (\ + 2[i)u xx exceeds k/1 2 , where I 
is the length of the stretched region. Hence, in graphene, 
stresses such that u xx ~ 10 -2 , will lead to wrinkles when 
I > lnm. The resulting deformation is periodic, with a 
wavelength, l w , and amplitude, A, in the direction per- 
pendicular to the applied strain given bv[4fj|: 



I h 



A = 



v/2 



1/4 



Aw 



,1/2 
1/2 



"11 



(6) 



Typical wrinkles, calculated using eq. ^ for a sus- 
pended region of length / = 200nm, and an applied strain 
u = 2 x 10~ 2 are shown in Fig. [5] 

A suspended sheet with constant tension gives rise to a 
constant gauge field 2 7]. Electron scattering takes place 
at the boundaries between the region with and without 
tension, and the the carriers propagate freely within the 
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suspended region. Wrinkles modulate the strains in the 
direction of the applied tension and also in the direction 
perpendicular to it. As a result, an effective magnetic 
field is induced throughout the suspended region. This 
field is also shown in Fig. [5l Its average value scales as: 

- c(3A 2 <f> cPu$o 



Electronic transport in suspended graphene samples is 
expected to be close to the ballistic limit [541] . Hence, even 
the low fields induced by wrinkles may lead to observable 
effects. 

Effects of effective gauge fields on Landau levels. The 
gauge fields induced by stresses modify the Landau levels 
induced by a magnetic field, and give to them an index 
dependent width, see Sec. VIII in The energy of 

the n-th Landau level scales as e„ oc \fn~B where B is 
the magnetic field. Stresses modify the effective field so 
that, to a first approximation the width of the Landau 
level becomes proportional to s/n. This effect is par- 
tially compensated by the averaging of the effective field 
within the spatial extent of the Landau level, given by 
a scale oc y/nlB, where Ib = \/&o/(nB) is the magnetic 
length, and $o is the quantum unit of flux. The narrow- 
est Landau level is that with n = 0, in agreement with 
experiments [55j|. 

Conclusions. We have analyzed possible deformations 
of graphene layers, and the effective gauge fields that 
they induce. The shape of the layer can be fixed by the 
corrugated substrate to which it is pinned. If the defor- 
mations induce substantial stresses, wrinkles can be pro- 
duced spontaneously. A similar elastic instability leads 
to wrinkles in suspended graphene samples. 
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